Extended RPA with ground-state correlations 
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We propose a time-independent method for finding a correlated ground state of an extended time- 
dependent Hartree-Fock theory, known as the time-dependent density- matrix theory (TDDM) . The 
correlated ground state is used to formulate the small amplitude limit of TDDM (STDDM) which 
is a version of extended RPA theories with ground-state correlations. To demonstrate the feasibility 
of the method, we calculate the ground state of ^^O and study the first 2^ state and its two-phonon 
states using STDDM. 



I. INTRODUCTION 



The study of unstable nuclei is a subject of current experimental and theoretical interests. Self-consistent theories 
such as the Hartree-Fock Bogoliubov theory (HFB) and the quasi-particle random-phase approximation (QRPA), 
which have extensively been used for stable nuclei, have also been applied to unstable nuclei, and the importance 
\^ of ground-state correlations has been demonstrated Q, |^ Q . Introducing a pairing field, HFB and QRPA deal with 
pairing correlations in the framework of a mean field theory. In contrast to HFB and QRPA, the time-dependent 
density-matrix theory (TDDM) , which is one of extended time-dependent Hartree-Fock theories, deals with ground- 
state correlations as genuine two-body correlations. TDDM has been applied to giant resonances in stable nuclei 0,0 
and also to low- lying collective states in unstable nuclei 0, 1^- The small amplitude limit of TDDM (STDDM) 
. 0, which is a time-independent version of TDDM, has also been used to calculate low-lying states in an oxygen 
T— I ' isotope ^3 ■ The importance of ground-state correlations has also been demonstrated by these TDDM and STDDM 
^-H calculations. The correlated ground state used in these nuclear structure calculations, however, is an approximate 
one which is obtained using a time-dependent method: The initial Hartree-Fock (HF) ground state is evolved in 
time using the equations of motion in TDDM and a time-dependent residual interaction whose strength gradually 
approaches its intended value with a time constant r. In the case of a solvable model where we can take r quite large, 
the ground state obtained has been found practically stationary and close to the exact one|llj|. However, in realistic 
cases where it is difficult to take sufficiently large r, the mixing of excited states, which causes spurious oscillations 
^ ' of some ground-state quantities, is unavoidable though it is small 12j. Therefore, it is anticipated to develop another 
^ , method for finding a correlated ground state of the TDDM equations. In this paper we propose a time-independent 
approach based on Newton's gradient method and demonstrate its feasibility by calculating the ground state of ^^O 
for which we have previously performed time-dependent calculations (SiiSj- ^^O is one of neutron-rich nuclei which 
attracts recent experimental and theoretical interests, and is quite suitable for our present study: Although it is an 
^ I open shell nucleus, the HF assumption which we use to obtain a starting ground state in the gradient method is valid 
as first-order approximation, and the omission of the proton degrees of freedom may be allowed in this explorative 
study because of the proton shell closure. The obtained ground state is used to construct the hamiltonian matrix of 
STDDM, and the first 2+ state and its two phonon states in ^^O are calculated. The paper is organized as follows. 
In sect. 2 the time-independent method for obtaining a correlated ground state is presented. In sect. 3 the formalism 
of STDDM is given. The results of numerical calculations for the ground-state, the first 2+ state, and its two phonon 
states are shown in sect. 4 and sect. 5 is devoted to a summary. 

II. METHOD FOR FINDING A CORRELATED GROUND STATE 

The ground state |$o) in TDDM is constructed so that 

Fiiaa') = ($o|[4aa',i?]|*o) = 0, (1) 



F2iaia2a2a[) = {'i>o\[a+^a+^aa'^ac,'^, H]\<i>o) = 0, (2) 

are satisfied, where H is the total hamiltonian consisting of the kinetic energy term and a two-body interaction, 
and [ ] stands for the commutation relation. In other words, the occupation matrix n^^, = {^o\a'^,aa\^o) and the 
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two-body correlation matrix C°^^^^, ^, = {^Q\a'^, a^, aa^o-ax l^'o) ~-^("-qiq' '^qjo' '^'^sre A is the antisymmetrization 
operator, are determined so that Eqs.lQJ and ^ are satisfied. The expressions for Eqs.lQJ and l(2Jl have already been 
given in Ref.|3 but are shown again in Appendix A. The single-particle wavefunction ipa is chosen to be an eigenstate 
of the mean field hamiltonian hQ{pQ): 

/i0(P0)^a(l) = -^7r^V'a(l) + / d2i; (1 , 2) [po (2, 2) Va ( 1 ) " PO ( 1 , 2)^Aa (2)] = Ea^a ( 1) , (3) 



2rn 

where the numbers denote space, spin, and isospin coordinates, and the one-body density matrix po is given as 

po(ll') - (4) 

aa' 

Attempts have been made to find a solution of Eqs.(^ and ^ However, this is not evident partly because 

Eqs.l^ and ^ are not in the form of an eigenvalue problem. The time-dependent method has been developed and 
tested for a solvable model JJj and realistic nuclei J, 8, 12] as mentioned above. In the following we propose a 
time-independent approach using iterative Newton's gradient method. We start from the HE ground state where 
^aa' = ^aa' (0) for occupicd (unoccupicd) single-particle states and C^^^^^' a' ^ 0. Then we iterate 



n"(iV + 1) \ _ / n"(Af) \ / 5Fi/5n^ 5Fi/SC° \ ^ ( Fi{N) 
C^iN + 1)) - [ C0(7V) ) "1^ 5F2/Sn^ SF^/SC' ) \ F2{N) 

(^-"W cy(F,{N)\ 

until convergence is achieved, where the matrix elements a, b, c, and d, which also depend on the iteration step N, 
are equivalent to those appearing in the hamiltonian matrix of STDDM. They are given in Ref.|lOl| and also shown 
in Appendix B. We have to introduce a small parameter a to control the convergence process. We have tested this 
iterative method for a solvable model and found that the obtained result is equivalent to the solution which had 
been obtained using the time dependent approach [Tl|. 



III. SMALL AMPLITUDE LIMIT OF TDDM 



TDDM gives the time-evolution of the one-body density-matrix p and the correlated part C2 of a two-body density- 
matrix P.ll5j. and STDDM has been formulated by linearizing the equations of motion for p and C2 9]. The equations 
of STDDM for the one-body amplitude Xaa'{p) and the two-body amplitude Aq^q^q'^q^ (/i) can be written in matrix 
form |1Q|| 

Eq.® can also be obtained from the following equations: 

{^o\K,a^,Hm - u^{^o\a+a^\<^), (7) 
(^o|[a+,a+ a„2a„i,i?]|$) = ujf,{'^^\a+a+ao,^aa,\<^), (8) 

where |$) is the wavefunction for an excited state with excitation energy w^. Linearizing Eqs.|(7I) and © with respect 
to Xaa' = (^ola^'Oal^) and Xa^a2a'^a'^ = (^0 1 o^^^ I 7 and using n°„, = ($o|aiaa|$o) and C^^^^„'^„^ = 
($o|aQ' aQ2'^ai I'^'o) — -^("-qiq' '^aaaOj Can obtain Eq.©. The fact that the linearization of Eqs.JTj) and © 
gives Eq.® might explain why the matrices a, b, c, and d are given by the variation of Fi and F2. The hamiltonian 
matrix of Eq.© is not hermitian as easily be understood from its explicit form (see Appendix B). In the case of 
a non-hermitian hamiltonian matrix, the ortho-normal and completeness relations are given not by the hermitian 
conjugate of \p) ~ {x, X) but by the left-hand-side eigenvectors \fi) = (x, X) which satisfy 



(9) 
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as explained in Ref.l^. When the ground-state |$o) is assumed to be the HF one and only the particle (p) - hole 
(h) and 2p - 2h amplitudes (and their complex conjugates) are taken in Eq.®, STDDM reduces to the second RPA 
(SRPA) [illlllll. 

The strength function is defined as 

SiE) = \{^,mo)\'S{E-E,), (10) 

E^>0 

where Iv^o) is the ground-state, l^*;^} is an excited state with an excitation energy E^, and Q an excitation operator. 
The strength function in STDDM for a one-body operator Qi is given as 



S{E) =- Um{ [[^(a|Qi|«'K'a(A*)) [YmiW^'M 



1 



E -uj^ + ir/2 



where an artificial width F is put to obtain a smooth distribution for S{E) and i^^, (/i) is defined as 

iL'(M) - • ^^')iAA'(A*) + ^i2(aa' : \i\2\\\'^)X^,x,x>^x>^{y). (12) 

AA' Aj^A2A-j^A?j 

Here and T12 are defined as 

Sii{aa' : W) = {<i>o\[a+ao,,a+ax']\'^o), (13) 
ri2(Q!a' : AiAaAiA'a) = {^o\[al^,aa,: a{^a{^ax'^ax'^ ■]\^o), (14) 

where : : means that : aX^aX^ax'^ax'^ := al^aX^ax'^ax'^ - A[aX^ax'^{<^o\aX^ax'^\^Q) + aX^_ax'^(<^o\aX^ax'^\<^o))- Similarly, 
the strength function in STDDM for a two-body excitation operator Q2 is given as 



-Re(wj,)>0 yc(iQ2Ct'j 



S{E) ^- -Im{ 2^ [ 2^ (aia2|Q2|a'ia2)^aiQ^QiQ2(M) 



Y if^lP2\Q2\m)Xj,',p',p,f3M 



1 



E-uj^ + iT/2 



Y ("ia2|(92|aia2)^a;a^aia2(A') 

Y (/3i/?2|Q2|/3i/3^)l^;,^,,fe(A*) I ^ , (15) 

where ^^^„,„^„^(m) is given by 

^aia2aia^(A*) = ^1 ^21 (aia2a'ia2 : ^A')iAA' (/i) 
AA' 

+ Yj '5'22(aia2aia2 : AiA2A'iA2)XaiA2a;a^(m)- (16) 

Ai A2Aj^A2 

Here Ti\ and ^22 are defined as 

?2i(aia2aia2 : ^^') = a^^a^.^aa-^aa^ ■.,a\ax'\\^{y) ■ (17) 

5'22(aia2aia2 : AiA2AiA2) = (^oli: a+,^a+ Ocj^aai :, : a+ a+ OA^flA; tH^o)- (18) 

'S'li, 5*22: 2^12 and T21 are written in terms of v^^^i and C'^^cC^', and are given explicitly in Ref. 19]. The strength 
functions in STDDM are not guaranteed to be positive definite, as is easily understood from Eas. (|ll|) and (|15|) . In 
practical applications this has not caused problems. 
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IV. NUMERICAL SOLUTIONS 
A. Correlated ground state 

In this subsection we present how the correlated ground state of ^^O is calculated and discuss some properties of 
the obtained ground state. To prepare the starting ground state for Eq.jSJ, we perform a static HF calculation as the 
first step. The Skyrme III (SKIII) is used as the effective interaction to generate a mean field. It has often been used 
as one of the standard parameterizations of the Skyrme force in nuclear structure calculations even for very neutron 
rich nuclei 0, HO]- We assume in the HF calculation that 1^5/2 is the last fully occupied neutron orbit of ^^O. Since 
Eq.(|3Jl involves matrix inversion at each iteration step, it is difficult to use a large number of single-particle states 
to solve Eq.lsj. In this explorative study we only use the two neutron orbits, 1^5/2 and 2si/2j to evolve 'nP^^i and 
^apa'i3'- Configurations consisting of these two orbits are considered to be the major components of the correlated 
ground state and the first 2+ state. The single-particle wavefunctions are confined to a sphere with radius 12 fm. 
The mesh size used is 0.1 fm. At each iteration step of Eq.©, Eq.lO is solved. This means that the single-particle 
states also evolve in a consistent way. In a self-consistent calculation the residual interaction should be the same as 
that used to generate the mean field. However, we had found in previous studies for giant resonances [23| that SKIII 
induces negligible ground-state correlations when the single-particle space is significantly truncated. This is the case 
also in this study of ^^O as will be shown below. Therefore, we p resent the results using the following pairing-type 
residual interaction of the density-dependent 5 function form 22] which is known to induce significant ground-state 
correlations Q 

v{r - r') = wo(l - p{r)lpo)5^{r - r'), (19) 

where p(r) is the nuclear density. The parameters po and vq are set to be 0.16fm~'^ and —900 MeV fm'^, respectively. 
Similar values of po and Vq have been used in HFB calculations [2^ 0, in truncated single-particle space. The 
converged result somewhat depends on how the iteration process is chosen. We found that to achieve real convergence 
it is necessary to start with a small value of vq and gradually increase it: In the calculation shown below we start 
with Wo/20 and increase it by Wo/20 for each 100 iterations. The sum of the absolute values of the matrix elements of 
Fi and F2 is shown in Fig.l as a function of the number of iterations. The final value of | -f IF2 1) at 2100th step 

is 2.5 X lO^'^MeV. Figure 1 demonstrates that the iteration method Eq.(|3J) works quite well. The final single-particle 
energies of the 1^5/2 and 2si/2 are -6.9 MeV and -2.8 MeV, respectively. Their occupation probabilities are 0.98 and 
0.05, respectively. The correlation energy i?cor defined as 

Eco. = \ («/3|^^|a'/3')Ca'/3'a/3 (20) 

is -0.60 MeV. The total energy is decreased by 0.22 MeV from the HF value: The correlation energy is largely 
compensated by the increase of the HF energy. 

B. Low- lying 2+ states 

Using the correlated ground state shown above, we solve Eq.® and calculate the strength function for 2+ states 
according to Ea. ((TT|l . To be consistent with the calculation of vP^^, and C^^^,^,, we only use the neutron lc?5/2 and 
2si/2 orbits for Xaa' and Xapa'fS' ■ The eigenvalues of some 2"*" states become imaginary because the hamiltonian matrix 
of Eq.(|SJl is not hermitian. However, their imaginary parts are quite small (less than 0.05MeV). Some 2+ states have 
also negative quadrupole strengths because positivity of S{E) is not guaranteed. However, the negative contributions 
are so small that S{E) becomes positive in the entire energy region when it is smoothed with F = 0.5MeV. The 
obtained result for Qi = r^l20 in STDDM (solid line) is shown in Fig. 2, where the strength functions in RPA (dot- 
dashed line) and SRPA (thin dotted line) is also presented for comparison. The 2+ state calculated in STDDM is 
energetically shifted upward and becomes significantly more collective as compared with that in RPA. The increase in 
the excitation energy is due to the lowering of the ground state, which is realized by the increase in the unperturbed 
p-h energy through the coupfing to C°,^,^^. We will discuss this point in more detail below. The enhancement of 
the collectivity of the 2+ state is due to the mixing of two-body configurations. These properties of the first 2+ state 
under the influence of ground-state correlations are similar to those obtained from QRPA calculations 0, Q . The 
reduced transition probability obtained is 68 fm"*. If the neutron effective charge is assumed to be 0.5e, the value 
of B{E2 : 0+ 2^) becomes 17 e^fm"*, which might be comparable with the experimental value of 21 ± 8 e^fm** 
[2^ . The 2+ state in SRPA is simply shifted downwards due to the coupling to 2p-2h configurations located around 
8 MeV. 
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FIG. 1: Sum of the absolute values of the matrix elements of Fi and F2 as a function of the number of iterations. Each spike 
corresponds to an increase of the interaction strength by vo/20 (see text). 




FIG. 2: Strength distributions of the neutron quadrupole modes in O calculated in STDDM (solid line), SRPA (thin dotted 
line) and RPA (dot-dashed line). The strength functions are smoothed with F = 0.5 MeV. 



C. Two phonon states 

The strength functions for the two-body operators Q2 = [r'^Y2<S^r'^Y2]Q with J =0, 2, and 4 are also calculated using 
Ea. (|15|l and the obtained results are shown in Fig. 3. The transition strength decreases with increasing J, whereas it 
would be independent of J if a pure phonon picture were valid. The ratios of the transition strengths of the 2+ and 
4+ states to that of the 0+ state are 0.55 and 0.30, respectively. These ratios should be compared with 4/5 and 1/3, 
respectively, which are obtained assuming unperturbed two 2si/2 particle - two 1^5/2 hole configurations. Although 
the excitation energies of the 2+and 4+ states are about twice the excitation of the 1 phonon state, the 0+ state is 
located at quite low energy. This originates in the nature of the residual interaction of the (5-function form (Ea. H19|) ) 
which strongly favors a nucleon pair with J'^ = 0+. The large spacing between the 0+ state and the 2+ state is 
similar to that obtained from the shell model calculation The strength distribution of the 2+ state has some 

small components in the low energy region (below 5 MeV) . This problem will be discussed in the next subsection. The 
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FIG. 3: Strength distributions of tlie two-phonon states of the neutron quadrupole modes in O calculated in STDDM. The 
solid, dotted and dot-dashed lines depict the results for 0"*", 2^ and 4^ states, respectively. The strength functions are smoothed 
with r = 0.5 MeV. 

transition probability between the one-phonon state and one of the two-phonon states is not well defined in STDDM 
but may be calculated as 



The reduced transition probabilities from the 0+, 2+, and 4+ states to the 1 phonon state are 16, 14, and 12 fm"*, 
respectively, which are close to the value 68/5«14 fm"* for the transition from the firs 2+ state to the ground state. If 
a pure phonon picture were valid, these values would be independent of the angular momenta of the 2 phonon states. 

In the following we compare STDDM with SRPA. The strength functions calculated in SRPA are shown in Fig. 4. 
The transition strengths in SRPA are much smaller than those in STDDM. Note that Fig. 4 is drawn in the same 
scale as Fig. 3. It had been pointed out _28 ] that the Xphph amplitude is important to reproduce the collectivity of 
low-lying two phonon states. To investigate this point, we performed a calculation for the 0+ state using a modified 
SRPA (mSRPA) which includes Xphph- mSRPA can be obtained from Eq.(6) by evaluating the hamiltonian matrix 
using the HF ground state. The result is shown in Fig. 5. The Xphph amplitude significantly enhances the collectivity 
of the 0"*" state. However, the transition strength in mSRPA is still smaller that that in STDDM. In Fig. 5 we also 
show the result of a further modified SRPA which includes all two-body amplitudes except for the 3p-lh, lh-3p, 3h-lp 
and lp-3h amplitudes. The 3p-lh and 3h-lp amplitudes do not increase the transition strength, although they change 
somewhat the form of the strength function. Now the transition strength becomes 91% of the strength in STDDM. 
Thus all two-body amplitudes except for the 3p-lh and 3h-lp types seem to be important to describe the collectivity 
of the low-lying 2 phonon states. This is in contrast to two-phonon states of giant resonances where the transition 
strengths are exhausted by the 2p-2h, 2h-2p and Iplh-lplh amplitudes 29]. The fact that the 0"^ state in STDDM 
is not much lowered as compared with the modified SRPA results is due to ground-state correlations. 

We also performed a self-consistent STDDM calculation using SKIII as a residual interaction. For simplicity, we 
neglected the spin-orbit force. Since there is a strong cancellation between the momentum dependent and independent 
terms, SKIII induces quite weak ground-state correlations in the truncated single-particle space considered in this 
work: The occupation probabilities of the 1^5/2 and 2si/2 are 0.9999 and 0.0003, respectively. The strength functions 
for the two-phonon states are shown in Fig. 6. The ratios of the g.s.^ 2+ and 4+ transition strengths to the g.s.^ O'*" 
one are 0.77 and 0.33, respectively, which are close to 4/5 and 1/3 obtained assuming the unperturbed configurations 
as discussed above. Since SKIII is effectively a very weak interaction in the truncated single-particle space, there is 
no significant difference between the STDDM and SRPA results. 



{<^o\[0^,Qi]\f/){^o\[0^,,Qi]\fi) 



(21) 



where is defined as 




(22) 
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FIG. 4: Strength distributions of the two-phonon states of the neutron quadrupole mode in O calculated in SRPA. The solid, 
dotted and dot-dashed lines depict the results for 0"*", 2"*" and 4^ states, respectively. The strength functions are smoothed 
with r = 0.5 MeV. 
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FIG. 5: Strength distributions of the 0"*" two-phonon state of the neutron quadrupole modes in ^^O calculated in STDDM (solid 
line), SRPA (dot-dashed line) and mSRPA (dotted line) which includes Xphph type amplitudes. The thin solid line indicates 
the result calculated using further modified SRPA where the Xpppp and Xhhhh type amplitudes are added to mSRPA. The 
strength functions are smoothed with F = 0.5 MeV. 

D. Incoherent states 

Now we discuss incoherent states which are seen in the low energy region of the strength function of the 2 phonon 
state with J'^ = 2'^ (Fig. 13 . These incoherent states have large components of the 3p-lh and lp-3h configurations 
which have the same unperturbed energy as the Ip-lh configurations. Since the 3p-lh, 3h-lp amplitudes and their 
complex conjugates are used as dynamical ampHtudes in Eq.@, these incoherent states naturally appear as eigenstates. 
However, the important physical role played by these amplitudes seems to be a rather kinematical one. To clarify this 
point, we use a different expression for Eq.® When the eigenvector (x, X) in STDDM is transformed to (y, Y) 
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FIG. 6: Strength distributions of the two-phonon states of the neutron quadrupole modes in O calculated in STDDM using 
SKIII. The solid, dotted and dot-dashed lines depict the results for 0^, 2"*" and 4"*" states, respectively. The result in SRPA for 
the 0^ state is also shown with the thin solid line. The strength functions are smoothed with F = 0.5 MeV. 



as 



Eq.® becomes 



X \ _ ( Sii Ti2 \ ( y 
X [T21 S22 \Y 



aSii + CT21 aTi2 + CS22 \ ( v \ _ f Sn Ti2 \ f y 

bSii + dT2l bT,2 + dS22 [Y '^''[T21 S22 [y 



(23) 



(24) 



It has been pointed out that this form of STDDM is a reasonable approximation for an extended RPA with 
hcrmiticity |19|. The kinematical effects of ground-state correlations arc included in the hamiltonian matrix of 
Ea. (|24|l . For example, the increase of the unperturbed energy of 1 phonon states due to ground-state correlations 
is given by CT21 in aSu + CT21, whereas aSu describes the renormalization of the RPA matrix due to the change 
in occupation factors. Omitting the 3p-lh and 3h-lp amplitudes and their complex conjugates in Y , we can avoid 
dynamical contributions of these amplitudes. This modified STDDM is referred to as mSTDDM. The obtained 
results for the 1 phonon state and the 2 phonon state with J'^ = 2+ are shown with thin dotted lines in Figs [7| and 
IHl respectively. The coherent states are little affected by the omission of the 3p-lh and 3h-lp amplitudes and their 
complex conjugates, and the incoherent states seen in the 2 phonon state can be eliminated. Thus it would be better 
to use Eg. 1)24(1 without the 3p-lh and 3h-lp type amplitudes instead of using Eq.© with full two-body amplitudes 
when we are interested in the calculation of the strength functions. However, we also found that the omission of the 
the 3p-lh and 3h-lp type amplitudes affects the transition probability between excited states when it is evaluated 
using Eq.lOJ. 

Finally we discuss incoherent states associated with one-body and two-body amplitudes which have zero unperturbed 
energy: Xpp, Xhh, Xphph, Xpppp, and Xhhhh are such amplitudes. Some incoherent solutions of Eq.@ or Ea. (|24|l have 
finite energies but most of them stay at zero or nearly zero energies. The transition strengths to these states are so 
small that they are invisible in the strength functions shown above. If these amplitudes are neglected and only Xph, 
Xhp, Xpphh, and Xhhpp are taken in Eg. 1)24(1 . these incoherent states disappear. However, a serious problem arises as 
to the collectivity of coherent states, especially of 2 phonon states, as explained above. Therefore, what we must do 
would be to interpret these incoherent states as unphysical ones and consider only coherent states. 
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FIG. 7: Strength distributions of the neutron quadrupole mode calculated in STDDM (solid line) and mSTDDM (thin dotted 
line). The result in RPA (dot-dashed line) is also shown for comparison. The strength functions are smoothed with F = 0.5 
MeV. 
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FIG. 8: Strength distributions of the 2+ two-phonon state calculated in STDDM (solid line) and mSTDDM (thin dotted 
line). The strength functions are smoothed with F — 0.5 MeV. 



E. Effects of three-body amplitudes 



The hamiltonian matrix 6T12 + ^5*22 for the two-body amplitudes in Ea. H24l) cannot represent all terms in the 
extended RPA of Ref . [l^ . The terms coming from 

^ e(Q;iQ;2a'ia2 : 71 727371 72 73)^32 (71 727371 7273 : ^I'^zAiAj) (25) 

7172737(7273 

are the missing terms. Here e is a matrix that would come into the hamiltonian matrix if we included a three-body 
amplitude in the evaluation of the left-hand side of Eq.®, and is given as 

e{aia2a[a2 : 7i72737i7273) = -{ail3\v\lil2)Sa^jJa'^Y^Sa'^^^ + (a273|i'l7i72)5ai73'5„;^j(5„^^^ 

+ (7l72b|ai73)'5ai7l'5a272'5ai73 " (7l72 ^ I^Ws) '5ai7l '^"272 '5a',7^ • (26) 
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FIG. 9: Strength distributions of the 0^ two-phonon state calculated in mSTDDM with Eg. 12511 (solid line) and without it 
(dotted line). The strength functions are smoothed with F = 0.5 MeV. 



is defined as 

732(7172737^7273 : AiA2AiA2) = ($o|[: a+ a+ a+, a.^^ a^, :, : a+ a+ aA^aA'J|$o)- (27) 

This eT^2 term contributes to the two-body sector of the hamiltonian matrix in Ea. (|24|l when the transformation 
of Eq.|(2SJ includes the three-body amplitudes. Kinematical effects of ground-state correlations such as the increase 
in unperturbed energies of 2 phonon states are expressed by some terms in Ea. H25|) . We performed a calculation in 
mSTDDM including Ea. H25() and found that the 1 phonon state is little affected: The excitation energy of the first 2+ 
state is unchanged and the transition strength to the first 2+ state is slightly increased (by 2.5%). The inclusion of 
the terms in Ea. H25|) somewhat affects the properties of the 2 phonon states as expected: The excitation energies are 
increased by 0.2 ^ 0.5 MeV and the transition strengths are decreased by 8 ~ 20 %, depending on J of the 2 phonon 
states. As an example the strength functions for the 0+ state calculated in mSTDDM with the terms in Ea. H25|) 
(solid line) and without them (dotted line) are shown in Fig. Thus the terms in Ea. H25|) need to be included in 
quantitative study of 2 phonon states. 



V. SUMMARY 



We proposed a time-independent method for obtaining a correlated ground-state of the time-dependent density- 
matrix theory (TDDM). The method was applied to obtain the correlated ground state of ^^O. The eigenstates of the 
small amplitude limit of TDDM (STDDM) were calculated for the first 2+ state in ^^O using the correlated ground 
state. It is found that STDDM properly deals with the effects of ground-state correlations on the low-lying 2+ state 
and that the non-hermiticity of STDDM is quite moderate one: The eigenvalues have quite small imaginary parts 
and the strength function is practically positive definite although it is not guaranteed in its non-hermitian form. The 
2 phonon states of the 2+ state were also studied. It was found that the 2 phonon state with J"^ = O"*" appears at 
very low excitation energy. This originates in the nature of the zero range force used. The results obtained using 
the Skyrme III force as a residual interaction were also presented. It was found that SKIII acts as a weak residual 
interaction in the very truncated single-particle space considered in this study. The physical roles played by the 3 
particle -1 hole and 3 hole - 1 particle type amplitudes in STDDM were discussed and a method for eliminating the 
incoherent states associated with these amplitudes was presented. It was also pointed out that the self-energy terms 
for unperturbed 2 phonon configurations, which are missing in STDDM, can be included using the extended RPA 
formalism of Ref . ,10. .IQJ . 
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APPENDIX A 

When tpa is chosen to be an eigenstate of the mean field hamiltonian (Eq.Q), Eqs.lQJ and 10) become 

(ea'-eo)nL. - E (^MA.a'A3(«A3|«|AiA2) - C^3^,;,, (AiAs^la'Aa)) (Al) 
Ai A2 A3 

where 

■B°/3a'/3' = (^1^2|w|A3A4)A[('5aAi -<Ai)(^/3A2 -»^l^Aj'^A3a'»^A4/3' 

Al A2 A3 A4 

- n\nl^^{6x^a' - nl,^^,){Sx^l3' ~ nl^fi,)], (A3) 

Paf3a't3' = X! (^1 ^2 I I A3 A4) [((^aAi <5;3A2 " ^^qAi '^/J A2 ~ '^Q a//3A2 )C'a3 A4a'/3' 
Al A2A3 A4 

- {Sx,a'Sx^p> - 5A3a'"A4/3' " '^A3a' ^A4/3' )C°/3Ai A2 ] ' (A4) 

Haf3a'f3' = ^ ( Al A2 | W | A3 A4) A [<5q Ai ("-Aga' C'a4/3A2/3' ~ "-A3/3' C'A4/3A2a' ) 
Al A2A3 A4 

+ '^/3A2("-A4/3'^A3aAiQ' ~ "■A4Q'^A3qAi/3') 

- <^q'A3 ('^aAi C'a4/3A2/3' ^ ""/SAi C'a4qA2/3' ) 

- <^/3'A4(«-53A2^A3qAiq' ~ "■aA2'-^A3/3Aia')]- (^5) 

Here the subscript A indicates that the corresponding matrix is antisymmetrized. 

APPENDIX B 

The matrices in Eq.© arc shown below: 

a(aa' : AA') = {ea - ea')Sa\Sa'X' -Y^{{f3X'\v\a X)An^f3 ~ {aX'\v\/3X)An"0a')i (Bf) 

13 

b{aia2a[a2 : AA') = -Sa^xi^HSa^p - <2/3)"7ai'^"a^ + "qs/jI'^t^'i " n°c[)iSsa'^ - ^s^'JiX' f3\v\-iS)A 

+ E((A'a2|^;|/37)C°,„,„, + {X' P\v\a[j) aC^.^^^p - {X' P\v\a',j) aC^.^-^p} 

Pi 

+ '5a2A{E[(<5ai/3 " <i/3)«°o'i + </3(Vi " "7«'i)('^*"2 " "L^ ) (^'/3|^'l7'^) a] 

MS 

+ E[(^'«ilH/37>C^,oi^ + (A'/3k|aW)AC°^,„,^ - (A'/3|T;|aW)^C°^,„,^]} 

+ '5a;A'{5I[('^'5"2 ~ "Li)"Qi/3'^a27 + "^^('^"i/' " ^°i/3)('5a27 " ""27) (/^^^l A(5) 

+ 5][(/37k|Aa^)C0^„,^^ + («i/3b|A7)AC^„,^ - {a^P^X^) aCI^,^,J} 

- '5aiA'{5I[('^'5"i - "L;)"ai/3'^a27 + (-^ai/S " ^^ai/s) ('5a27 " "027) (/^^l^^l A(^) a] 

/37'5 

+ E[(/37lHAa;)C°^„,^^ + (ai/3|^;|A7)AC^„,0 - (a2/3|t;|A7)AC°^^„,^]} 
+ E[(«i^>I/3A)aC^„,„;„^ - (a2A>|/3A)^C°„^„,„, 

13 

- {f}X'\v\a'^X)ACl,^,^.^p + {pX'\v\a',X)ACl^^,^.J, (B2) 



c(aQ;' : A1A2A1A2) = {aX'2\v\XiX2)Sa'x'^ ~ {X'iX2\v\a'X2)da\i, 
d{aia2a[a'2 : AiA2A'iA2) = (e^i + - ea[ - ^a'J^aiXiSa2\2^a[\[Sa'^x'^ 

+ Sa'^x[Sa'^x'2 ^{Sail3Sa2f " '5a27"-ai/3 " ^ai pn^^^) {P'y\v\XiX2) 
fi't 

- Sa,xja2\2 X!(^"'i/3'^"^7 - '^a^7"-a;/3 " <^c;/3»^^^) ( A'l A'a | W |/37) 

+ ^o,2\2^a'2\'2^ii(^i^'i\'"\(^^i)Anl^,^ - {f3X[\v\a[Xi)An°^,0) 

P 

+ ^o,2\2^a[\[^{{(^i>'2\v\f3Xi)Anl^,^ - {PX'2\v\a'2Xi)An'i,p) 



+ Sa,XiSa'^x'2'^{{a2X[\v\pX2)Anl^,^ - (/3A'i |u|a'iA2)A<2/3)' 
where n°, = (^olaiaal'I'o) and C° , , = ($o|a^' a^' acfla, I'J'o) - , n° ,)• 



[1] E. Khan and N. Van Giai, Phys. Lett. B 472 (2000) 253. 

[2] E. Khan, N. Sandulescu, M. Grasso and N. Van Giai, Phys. Rev. C66 (2002) 024309. 

[3] M. Matsuo, Nucl. Phys. A696 (2001) 371. 

[4] M. Gong and M. Tohyama, Z. Phys. A335 (1990) 153. 

[5] F. De Blasio et al., Phys. Rev. Lett. 68 (1992) 1663. 

[6] M. Tohyama, Nucl. Phys. A657 (1999) 343. 

[7] M. Tohyama and A. S. Umar, Phys. Lett. B516 (2001) 415. 

[8] M. Tohyama and A. S. Umar, Phys. Lett. B549 (2002) 72. 

[9] M. Tohyama and M. Gong, Z. Phys. A332 (1989) 269. 

[10] M. Tohyama and P. Schuck, nucl-th/0305088 

[11] M. Tohyama, Prog. Theor. Phys. 92 (1994) 905. 

[12] M. Tohyama, Prog. Theor. Phys. 94 (1995) 147. 

[13] M. Tohyama, P. Schuck and S. J. Wang, Z. Phys. A339 (1991) 341. 

[14] H. J. Lipkin, N. Meshkov and A. J. Ghck, Nucl. Phys. 62 (1965) 188. 

[15] S. J. Wang and W. Gassing, Ann. Phys. 159 (1985) 328; W. Gassing and S. J. Wang, Z. Phys. 328 (1987) 423. 

[16] J. Sawicki, Phys. Rev. 126 (1962) 2231; J. Da Providencia, NucL Phys. 61 (1965) 87. 

[17] G. Yannouleas, Phys. Rev. C35 (1987) 1159. 

[18] S. Droidz, S. Nishizaki, J. Speth and J. Wambach, Phys. Rep. 197 (1990) 1. 

[19] M. Tohyama and P. Schuck, nucl-th/0305087 

[20] T. Otsuka, N. Fukunishi and H. Sagawa, Phys. Rev. Lett. 70 (1993) 1385. 

[21] M. Tohyama, Phys. Rev. C58 (1998) 2603. 

[22] R. R. Ghasman, Phys. Rev. C14 (1976) 1935. 

[23] J. Terasaki, H. Flocard, P. -H. Heenen, P. Bonche, NucL Phys. A 621 (1997) 706. 

[24] T. Duguet, P. Bonche and P.-H. Heenen, NucL Phys. A 679 (2001) 427. 

[25] M. Yamagami, K. Matsuyanagi and M. Matsuo, NucL Phys. A 693 (2001) 579. 

[26] P. G. Thirof et aL, Phys. Lett. B485 (2000) 16. 

[^27] B. A. Brown, Prog, in Part, and Nucl. Phys. 47 (2001) 517; http:/ /www.nscLmsu.edu/~br own/database.htm| 

[28] N. Kanesaki, T. Marumori, F. Sakata and K. Takada, Prog. Theor. Phys. 50 (1973) 867. 

[29] M. Tohyama, Phys. Rev. C64 (2001) 067304. 



